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A Correlation Between the Higgs Mass and Dark Matter 
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Depending on the value of the Higgs mass, the Standard Model acquires an unstable region at 
large Higgs field values due to RG running of couplings, which we evaluate at 2-loop. For currently 
favored values of the Higgs mass, this renders the electroweak vacuum only meta-stable with a long 
lifetime. We argue on statistical grounds that the Higgs field would be highly unlikely to begin in 
the small field meta-stable region in the early universe, and thus some new physics should enter 
in the energy range of order, or lower than, the instability scale to remove the large field unstable 
region. We assume that Peccei-Quinn (PQ) dynamics enters to solve the strong CP problem and, 
for a PQ-scale in this energy range, may also remove the unstable region. We allow the PQ-scale to 
scan and argue, again on statistical grounds, that its value in our universe should be of order the 
instability scale, rather than (significantly) lower. Since the Higgs mass determines the instability 
scale, which is argued to set the PQ-scale, and since the PQ-scale determines the axion properties, 
including its dark matter abundance, we are led to a correlation between the Higgs mass and the 
abundance of dark matter. We find the correlation to be in good agreement with current data. 



INTRODUCTION 



Recent LHC results are consistent with the predictions 
of the Standard Model, including the presence of a new 
boson that appears to be the Higgs particle with a mass 
rriH ^ 125 - 126 GeV [H, With the Higgs at this 
mass, the Standard Model is well behaved up to very 
high energies if we evolve its couplings under the renor- 
malization group (RG) equations. By no means does this 
imply that the Standard Model will be valid to these very 
high energies, and in fact there are good phenomenolog- 
ical reasons, such as dark matter, strong CP problem, 
baryogenesis, inflation, hierarchy problem, etc, to think 
it will be replaced by new physics at much lower energies, 
say C(TeV). But it is logically possible, albeit unlikely, 
that the Standard Model, or at least the Higgs sector, 
will persist to these very high energies and the explana- 
tion of these phenomena will be connected to physics at 
these high, or even higher, energy scales. 

So at what energy scale must the Standard Model 
breakdown? Obviously new physics must enter by the 
Planck scale Mpi where quantum gravity requires the in- 
troduction of new degrees of freedom. However, the RG 
running of the Higgs self-coupling A can dictate the need 
for new physics at lower energies, depending on the start- 
ing value of A. The Higgs mass is related to the self- 
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coupling by mH ~ v2At;£;n/, where the Higgs VEV is 
Vew ~ 246 GeV. For moderate to high values of the Higgs 
mass, the initial value of A, defined at energies of order 
the electroweak scale, is large enough that it never passes 
through zero upon RG evolution. On the other hand, for 
small enough values of the Higgs mass, the self-coupling 
A passes through zero at a sub-Planckian energy, which 
we denote E* , primarily due to the negative contribution 
to the beta function from the top quark, acquiring an un- 
stable region at large field values [s, [ij . The latter occurs 
for a light Higgs as has been observed. One finds that this 
renders the electroweak vacuum only meta-stable with a 
long lifetime. However, we will argue in this paper that it 
is highly unlikely for the Higgs field in the early universe 
to begin in the meta-stable region as that would require 
relatively small field values as initial conditions. Instead 
it would be much more likely to begin at larger field val- 
ues, placing it in the unstable region. Hence, the energy 
scale E* sets the maximum energy scale for new physics 
beyond the Standard Model to enter. 

There are many possible choices for the new physics. 
One appealing possibility is supersymmetry, which alters 
the running of the Higgs self-coupling due to the presence 
of many new degrees of freedom, likely entering at much 
lower energies, conceivably ©(TeV), or so. In addition to 
possibly stabilizing the Higgs potential, supersymmetry 
can also alleviate the hierarchy problem, improve unifi- 
cation of gauge couplings, and fit beautifully into fun- 
damental physics such as string theory. So it is quite 
appealing from several perspectives. It is conceivable, 
however, that even if supersymmetry exists in nature, it 
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FIG. 1: Dark matter density ( (in units of baryon density) 
as a function of the Higgs mass niH is given by the blue 
curves. The solid-blue curve is for the central value of the 
top mass mt = 173.1 GeV. The dashed-blue curves are for 
mt = 173.1 ± 0.7 GeV, with the upper value on the right and 
the lower value on the left. The red vertical lines indicate 
the measured Higgs mass range rriH = 125.7 ± 0.6 GeV from 
combining ATLAS and CMS data pj, |3| . The green horizon- 
tal lines indicate the measured dark matter density to baryon 
density ( — Qdm /^b range, where Udm = 0.229 ±0.015 and 
fis = 0.0458 ± 0.0016, from WMAP7 data 



is spontaneously broken at very high energies, and in such 
a scenario we would be forced to consider other possible 
means to stabilize the Higgs potential. 

One intriguing possibility that we examine in this pa- 
per is to utilize dynamics associated with the solution of 
the strong CP problem; the problem that the CP violat- 
ing term ~ d e^^'^"-^ F^^F^^ in the QCD Lagrangian is ex- 
perimentally constrained to have coefficient 16*1 < 10~^°, 
which is highly unnatural. The leading solution involves 
new Peccei-Quinn (PQ) dynamics involving a new 
complex scalar field (/) and a new global U{1) symmetry 
that is spontaneously broken at some energy scale FpQ. 
This leads to a new light scalar field known as the axion 
[1, 0]- Since it is bosonic, the (f) field adds a positive con- 
tribution to the effective A for the Higgs, potentially re- 
moving the unstable region, depending on the scale Fpq . 
This elegant mechanism to remove the unstable region 
was included in the very interesting Ref. [8|, where this 
and other mechanisms were discussed, and was a source 
of motivation for the present work (also related is 0). 



In the present paper, we would like to take this elegant 
mechanism for vacuum stability and push it forward in 
several respects. Firstly, as already mentioned, we will 
argue on statistical grounds why the meta-stable vacuum 
requires stabilization. Secondly, we will allow the PQ- 
scale to scan and argue, again on statistical grounds, why 
it should be of order the instability scale E* , rather than 
orders of magnitude lower. Finally, we will furnish a cor- 
relation between the Higgs mass and the axion dark mat- 
ter abundance, and use the latest LHC [H, 0] and cosmo- 
logical data [10] to examine the validity of this proposal. 
The outcome of this series of arguments and computation 
is presented in Fig. [1] which is the primary result of this 
work. 

The outline of our paper is as follows: In Section HIl we 
examine the running of the Standard Model couplings at 
2-loop order. In Section lTTll we examine the meta-stability 
of the Standard Model vacuum and argue that it is statis- 
tically unfavorable for the Higgs to begin in this region. 
In Section IIVI we include Peccei-Quinn dynamics to re- 
move the Higgs instability and argue why Fpg ^ E* . 
In Section |V] we relate the PQ-scale to the axion dark 
matter abundance, which furnishes a correlation between 
the Higgs mass and the abundance of dark matter. Fi- 
nally, in Section I VII we compare the correlation to data 
and discuss our results. 



II. STANDARD MODEL RG EVOLUTION 

We begin with a reminder of the structure of the Higgs 
sector of the Standard Model. The Higgs field is a com- 
plex doublet H with Lagrangian 



C = DH^'DH 



(1) 



In the unitary gauge we expand around the VEV as 
H = {0,vew + h)/^/2, where in our convention vew ~ 
246 GeV. The associated Higgs mass is mn — V2Xvew 
in terms of the starting value of A, normally defined 
around the Z boson mass. At higher energies, the self- 
coupling A undergoes RG evolution due to vacuum fluc- 
tuations from self interaction, fermion interactions, and 
gauge interactions. Defining dX/dt — I3\ with t — InE/^ 
the associated 1-loop beta function (suppressing external 
leg corrections for now) is 



(47r)2 



24A^ 



6yt + l{2g'+{g' + g'')') 



(2) 



where the only fermion Yukawa coupling we track is that 
of the top quark yt since it is by far the largest. For 
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sufficiently large Higgs mass, the positive self interaction 
term ~ is large enough to keep the beta function 
positive, or only slightly negative, to avoid A running 
negative at sub-Planckian energies. For sufficiently small 
Higgs mass, the negative top quark contribution ~ —yf 
can dominate and cause the beta function to go negative, 
in turn causing A to pass through zero at a sub-Planckian 
energy, which we denote E* . The top quark Yukawa cou- 
pling itself runs toward small values at high energies with 
1-loop beta function 
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which is quite sensitive to the value of the strong cou- 
pling gs ■ To compute the evolution of couplings and the 
quantity E* = E*{mH, J/t, • • •) accurately, we do the fol- 
lowing: (i) Starting with couplings defined at the Z mass, 
we perform proper pole matching and running up to the 
top mass, (ii) we include external leg corrections (and 
the associated wavefunction renormalization) , (iii) we si- 
multaneously solve the 5 beta function differential equa- 
tions for the 5 important couplings A, yt, g' , g, gs, and (iv) 
we include the full 2-loop beta functions for the Stan- 
dard Model; these are presented in the Appendix (see 
Refs. [H [H for more information). In our numerics, 
we use particular values of the couplings g\g,gs, derived 
from the best fit values 
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FIG. 2: Higgs self-coupling A as a function of energy, for 
different values of the Higgs mass from 2-loop RG evolu- 
tion. Lower curve is for mn = 116 GeV, middle curve is for 
TTiH = 126 GeV, and upper curve is for rriH = 130 GeV. All 
other Standard Model couplings have been fixed in this plot, 
including the top mass at mt = 173.1 GeV. 



III. META-STABILITY AND PROBABILITY 



^ TT^' sin^ Ow = 0.2311, a,(mz) = 0.1184. 

(4) 

In our final analysis, we will allow for three different val- 
ues of mt = V2yt VEW, namely the central value and 1- 
sigma variation mt = 173.1±0.7 GeV, and we will explore 
a range of mH = V^Xvew, with vew = 246.22 GeV. 

Performing the RG evolution leads to the energy de- 
pendent renormalized coupling X{E). A plot of X{E) 
is given in Fig. [2] for three Higgs mass values, namely 
mn = 116 GeV (lower curve), mn = 126 GeV (middle 
curve), and mn — 130 GeV (upper curve), with the top 
mass fixed to the central value mt = 173.1 GeV. This 
shows clearly that for the lighter Higgs masses that the 
coupling A passes through zero at a sub-Planckian energy 
scale E* and then remains negative. Furthermore, since 
the coupling only runs logarithmically slowly with energy, 
the value of E* can change by orders of magnitude if the 
starting value of the couplings changes by relatively small 
amounts. The domain E > E* involves a type of "attrac- 
tive force" with negative potential energy density, as we 
now examine in more detail. 



If we think of the field value h as being the typical 
energy pushed into a scattering process at energy E, then 
we can translate the RG evolution of the couplings into 
an effective potential. Using X{E) and replacing E h, 
we obtain the (RG improved) effective potential at high 
energies {h ^ vew) 

VM^\mG{trh\ (5) 

where the wavefunction renormalization factor G is given 
in terms of the anomalous dimension 7 by G(t) — 
exp(— /q 7(t')(it'), and we replace t — > ln/i//i. Hence 
for a Higgs mass in the range observed by the LHC, 
the effective potential Ktt goes negative at a field value 
h = E* that is several orders of magnitude below the 
Planck scale, as can be deduced from the behavior of 
A(£;) with mH = 126 GeV in Fig.H 

We could plot Kff (^) directly, however the factor of /i^ 
makes it vary by many orders of magnitude as we explore 
a large field range. Instead a schematic of the resulting 
potential will be more illuminating for the present discus- 
sion in order to highlight the important features, as given 
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FIG. 3: Schematic of the effective potential V^tt as a function 
of the Higgs field h. This is not drawn to scale; for a Higgs 
mass in the range indicated by LHC data, the heirarchy is 
Vew <S -E* <g Mpi, where each of these 3 energy scales is 
separated by several orders of magnitude. 

in Fig. [3l The plot is not drawn to scale; the 3 energy 
scales satisfy the hierarchy vew <^ E* <^ Mp, for a Higgs 
mass as indicated by LHC data rriH ^ 125 — 126 GeV. 
Note that the local maximum in the potential occurs at 
a field value that is necessarily very close to E* (only 
slightly smaller) and so we shall discuss these 2 field val- 
ues interchangeably. 

In this situation, the electroweak vacuum is only meta- 
stable. Its quantum mechanical tunneling rate can be es- 
timated by Euclideanizing the action and computing the 
associated bounce action ^o. This leads to the following 
probability of decaying in time Tu through a bubble of 
size R [13|] 

p^{Tu/Rre~s°. (6) 

The computation of the rate is rather involved, and we 
shall not pursue the details here. Suffice to say that for 
the central values of Higgs mass and top mass from LHC 
data, it is found that the lifetime of the electroweak vac- 
uum is longer than the present age of the universe fl3 . fl5} . 

It is conceivable that it is an acceptable situation for 
the electroweak vacuum to be meta-stable. However, here 
we would like to present an argument that such a situ- 
ation is statistically disfavorable. We imagine that in 
the very early universe, the Higgs field was randomly 



distributed in space. For instance, during cosmological 
infiation the Higgs field could have been frozen at some 
value as the universe rapidly expands (if high scale in- 
flation) until after inflation when the field will oscillate 
and its initial value could plausibly have been random 
and uniformly distributed. If this is the case, then what 
is the probability that the Higgs field began in the meta- 
stablc region h < E* , rather than the unstable region 
h> E*? The answer depends on the allowed domain the 
Higgs can explore. Here we estimate the allowed domain 
to be Planckian, i.e., < h < Mpi, but our argument 
only depends on the upper value being much larger than 
E* . Naively, this would lead to a probability ~ E* /Mpi, 
however we should recall that the Higgs is a complex 
doublet, composed of 4 real scalars, and each one would 
need to satisfy h < E* in the early universe to be in the 
meta-stable region. Hence, we estimate the probability 
as 

Prob (Higgs begins in meta-stable region) ^ ' 

(7) 

For instance, for niH ~ 125.5 GeV and mt = 173.1 GeV, 
we have E* ~ 10^^ GeV, leading to a probability ~ 
(10" GeV/lOi^CeV)"* = 10~^^, which indicates that the 
chance of randomly landing in the meta-stable region in 
the early universe is exceedingly unlikely. Instead it is 
far more likely to land in the unstable region indicated 
in Fig. [3l Here the effective potential is negative leading 
to a catastrophic runaway instability, perhaps to a new 
VEV that is close to Planckian. This would in turn lead 
to a plethora of problems for the formation of complex 
structures, etc, so we can safely assume such a regime is 
uninhabitable and irrelevant. This leads us to examine 
a scenario in which new physics enters and removes this 
problem. 

IV. PECCEI-QUINN DYNAMICS AND 
DISTRIBUTION 

One of the phenomenological reasons for new physics 
beyond the Standard Model is the fine tuning of the CP 
violating term in the QCD Lagrangian. The following 
dimension 4 operator is gauge invariant and Lorentz in- 
variant and should be included in the QCD Lagrangian 
with a dimensionless coefficient 9 

AC^^e'^-'-^F^^F^^. (8) 

From bounds on the electric dipole moment of the neu- 
tron, this term is experimentally constrained to satisfy 
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\9\ < 10""'^'', which requires extreme fine tuning. There 
appears to no statistical explanation of this fine tuning if 
it were purely random, since a small but moderate value 
of 9 would have very little consequences for the forma- 
tion of complex structures. Instead this requires a dy- 
namical explanation, which we take to be due to a new 
global symmetry, known as a Peccei-Quinn (PQ) sym- 
metry 5], involving a new heavy complex scalar field (f). 
This field undergoes spontaneous symmetry breaking at 
a scale FpQ , and in the resulting effective field theory, the 
quantity 9 is essentially promoted to the angular degree 
of freedom in^ = a light scalar field known as 

the axion 0] . The zero temperature Lagrangian for 
includes a symmetry breaking potential for p and a QCD 
instanton generated sinusoidal potential for 9 

C^\p\d,9f + \{d,pY 

-2A2sin2(0/2)-^(p^-F|Qf. (9) 

If we expand around the (/)-field's VEV at low energies, 
we see that the angular component of 0, the axion, is very 
light with a mass TOq = A^/Fpg (where A is of order the 
QCD-scale) and will be dynamically driven to zero, solv- 
ing the strong CP problem. Since we will require Fpq 
to be a very high energy scale (compared to say the elec- 
troweak scale), the radial mode of (p is very heavy, with 
mass FpQ. Hence at low energies, the radial mode is 
essentially irrelevant; it can be integrated out and, apart 
from a possible renormalization of the Standard Model 
couplings, can be ignored. However at energies approach- 
ing FpQ, the radial mode cannot be ignored. The field 
(j) couples to various Standard Model particles in any re- 
alization of the PQ symmetry, including the Higgs field 
through the interactions of the type ~ Wllcfi^cj). This 
causes an alteration in the effective coupling A at an en- 
ergy scale of order E ~ Fpq where the new field becomes 
dynamical. Since 4> is bosonic, it generally leads to a pos- 
itive increase in A, either through tree- level corrections or 
through loop corrections, leading to a threshold boost in 
Xa,t E ^ FpQ ; a point that was included nicely in Ref. [8] . 

This can lead to a reduction or removal of the unstable 
region depending on the scale Fpq relative to the in- 
stability scale E* , assuming 0(1) couplings between the 
Peccei-Quinn dynamics and the Higgs sector. We obvi- 
ously require Fpq to be in the range Fpq < E* in or- 
der for the new physics to prevent the effective potential 
V'off(^) from having a large negative regime (note that a 
small negative dip is statistically allowable for h, but not 
a large field dip). But since E* is very large, this leaves 
several orders of magnitude uncertainty in the value of 
FpQ. In other words, it would be sufficient for Fpq <^ E* 



in order to remove the unstable region. However, here we 
would like to present a statistical argument that 

FpQ ^ E* (10) 

is much more likely. There are indications that the PQ- 
scale may be associated with GUT or Planckian physics, 
and indeed typical realizations of the QCD-axion in string 
theory suggests that Fpq is much closer to the GUT or 
Planck scale Mpi [l^, rather than a more intermediate 
scale, such as ~ 10^^ GeV. In some landscape, we can 
imagine Fpq scanning over different values. For lack of 
more detailed knowledge, we can imagine that it scans on, 
say, a uniform distribution in the range < Fpq < Mp,. 
If this is the case, then Fpq will be as small as is required, 
but would not be significantly smaller as that would be 
even rarer in the landscape. By placing Fpq on a uniform 
distribution, the probability that it will be small enough 
to alleviate the instability is roughly 

E^ 

Prob(PQ-field alleviates instability) ~ — — , (11) 

Mpi 

where almost all of the phase space pushes Fpq ^ E* , 
rather than orders of magnitude lower. It is important 
to note that for E* ^ Afpi, as arises from the measured 
Standard Model's couplings, the probability in eq. (fTTj) 
is small but still much greater than the probability in 
eq. (I7|). Hence it is much more likely to have an atypically 
small PQ-scale and no constraint on the initial Higgs field, 
than an atypically small Higgs field and no constraint on 
the PQ-scale. We now examine the cosmological conse- 
quences of Fpq ^ E*. 

V. AXION DARK MATTER 

The light scalar axion particle is neutral, very stable, 
and acts as a form of dark matter. The computation 
of its abundance is non-trivial and has been studied in 
many papers, including Refs. [13, [3- The final result 
for the axion abundance is essentially controlled by the 
scale Fpq. Its value is normally measured by the quan- 
tity Qdm = RDM / Pent, where pdm is the energy density 
in axion dark matter and Pcrit is the so-called "critical 
density" of the universe defined through the Friedmann 
equation as Pcrit = iH'^ / {^i^G). Tracking the non-trivial 
temperature dependence of the axion potential and red- 
shifting to late times, leads to the following expression 

for ^DM 
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where the Hubble parameter is represented as H = 
100/ikm/s/Mpc and T is the CMB temperature. The 
coefficient x is an 0(0.1 — 1) fudge factor due to uncer- 
tainty in the detailed QCD effects that set the axion mass 
and its temperature dependence. In our numerics we have 
taken x — 0.5 as a representative value. It is quite possi- 
ble that the true value may be smaller than this, such as 
X ~ 0.15 as taken in Ref. [l^, but other effects, includ- 
ing contributions from string-axions, etc, can potentially 
push the true value to be larger p^]. Also, 9i is the ini- 
tial 6 angle in the early universe (which later redshifts 
towards zero, solving the strong CP problem). Here we 
take (6*1) = which comes from allowing 9i to be uni- 

formly distributed in the domain —TT<9i < tt and then 
spatial averaging. Another interesting possibility arises 
if inflation occurs after Peccei-Quinn symmetry break- 
ing, allowing 9i to be homogeneous and possibly small, 
as studied in Refs. [l^, HH. The latter scenario is subject 
to various constraints, including bounds on isocurvature 
fluctuations 
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and will not be our focus here. 
The quantity fl dm is slightly inconvenient for express- 
ing the main results for the following two reasons: (i) 
in a flat universe (as we are assuming) it is bounded to 
satisfy Q,dm < 1, which obscures the fact that a priori 
the dark matter abundance could be enormous, and (ii) 
it is manifestly time dependent (due to h and T), which 
requires some choice of physical time to compare differ- 
ent universes. To avoid these complications, we prefer to 
compute the dark matter density in units of the baryon 
density. Fixing the baryon to photon ratio at the mea- 
sured value, we have 



0.7 

h 



T 



2.725 K 



(13) 



with fls.o ~ 0.046 from observation. From this we define 
the (unbounded and time independent) measure of dark 
matter C as 





Pdm 




PB 




FpQ 




1012 GeV 



7/6 



(14) 



(15) 



Observations show that the dark matter density param- 
eter C, is non-zero in our universe, although its particu- 
lar particle properties (whether axion or WIMP, etc) are 
still unknown. The observational evidence for dark mat- 
ter comes from a range of sources, including CMB, lens- 
ing, galaxy rotation and clustering, structure formation, 
baryon-acoustic-oscillations, etc, and is very compelling, 
e.g., see Refs. (23l-l28|. and its abundance has been mea- 
sured quite accurately. Hence our prediction for the value 



of C (coming from setting Fpq ^ E* with E* determined 
by mn) can be compared to observation; see Fig. [1] 



VI. RESULTS AND DISCUSSION 

A. Comparison with Data 

Let us summarize our argument: Holding other pa- 
rameters fixed, the Higgs mass niH determines the insta- 
bility scale E* ^ which we evaluate at 2-loop. We have 
argued on statistical grounds in Section Hill whv the scale 
of new physics should not be larger than E* and in 
Section IIVI why the scale of new physics should not be 
(significantly) smaller than E* , leading to Fpq ^ E* . 
Since Fpq determines the dark matter abundance ^ in 
eq. (|15p . this establishes a correlation between niH and 
^. The result was displayed earlier in the paper in Fig. [TJ 
The solid-blue curve is for the central value of the top 
mass mt = 173.1 GeV, and the dashed-blue curves are 
for mt = 173.1 ± 0.7 GeV. We compare this prediction to 
the latest LHC and cosmological data. Firstly, we have 
takenthe ATLAS value mjy = 126.0±0.4(stat)±0.4(syst) 
l and the CMS value mn ^ 125.3 ± 0.4(stat) ±0.5(syst) 
2j, and produced our own combined value of mn — 
125.7 ± 0.6 GeV, which is indicated by the red verti- 
cal lines. Secondly, we have taken the WMAP7 data, 
plus other observations, for the dark matter abundance 
^DM = 0.229 ± 0.015 and the baryon abundance ^b — 
0.0458 ± 0.0016 JUj and combined to obtain C, which is 
indicated by the green horizontal lines. The predicted 
correlation between the Higgs mass niH and the dark 
matter abundance C in Fig. [T] displays good agreement 
with current data. 



B. Precision and Uncertainties 

Improved accuracy in testing this scenario comes in 
several experimental directions. This includes measuring 
the Higgs mass mn to better precision, as well as the top 
mass rrit and the strong coupling as (which we set to the 
central value as — 0.1184), while the current accuracy in 
C is quite good. A theoretical uncertainty surrounds the 
specific choice of Fpq relative to E* . Here we have taken 
Fpq ^ E* , due to a statistical argument that allowed the 
scale Fpq to scan, leading to the conclusion that it should 
be as small as required, but no smaller - an argument 
that is similar to the argument for the magnitude of the 
cosmological constant [23| . One might argue that a factor 
of a few smaller may be required to properly alleviate the 
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instability ^ , which would lead to a slight lowering of the 
blue curves in Fig.[TJ but a small negative dip is tolerable 
statistically, which makes Fpq ^ E* plausible. 

Related to this uncertainty is the particular prior dis- 
tribution for FpQ^ which we assumed to be uniform. The 
expectation of a flat distribution is plausible for the cos- 
mological constant A if one allows both positive and neg- 
ative values, making A ^ not special. In the case of 
FpQ, it is necessarily positive, so Fpq ^ is arguably a 
special part of the distribution. This may render the true 
distribution non-uniform. However, as long as the distri- 
bution does not vanish in the Fpq — limit faster than 
(Fpq)'^, then our arguments go through. In other words, 
the probability of an atypically small Fpq and no con- 
straint on the initial Higgs field would still be larger than 
the probability of an atypically small Higgs field and no 
constraint on Fpq. Also, one may question whether the 
uniform distribution assumed for the initial values of each 
of the 4 components of the Higgs field is reasonable. Since 
we have a sufficiently limited understanding of the early 
universe, including a measure problem for inflation, any 
such assumptions could be called into question. However, 
since the meta-stable region occupies such a tiny fraction 
of the volume of fleld space, roughly ~ (10~^)'* — 10^^^ 
or so, an alteration in prior probabilities would need to 
be quite drastic to change the conclusions. 

C. Outlook 

An important test of this scenario involves unravelling 
the nature of dark matter directly. The QCD-axion is 
actively being searched for in a range of experiments, in- 
cluding ADMX [s^j, with no positive detection so far. 
But the regime of parameter space in which the axion can 
be the dark matter will be explored in coming years. If 
an axion is discovered, it will be important to unravel its 
particular properties including its coupling to other flelds. 
An explicit embedding of the discovered version (popu- 
lar models include KSVZ P,'3^ and DFSZ 'slQ) into 
the Higgs stability analysis would be important. Searches 
such as ADMX rely upon the axion being all or most of 
the dark matter, so a related verification would be the as- 
sociated lack of discovery of WIMPs, or other dark mat- 
ter candidates, in direct or indirect searches. Or at least 
these forms of dark matter should comprise a relatively 
small fraction of the total. 



The discovery of the Higgs boson at the LHC is a fi- 
nal confirmation of the Standard Model. This leaves the 
scale at which the theory breaks down unclear. Here we 
have investigated the possibility that the theory, or at 
least the Higgs sector, remains intact until the scale at 
which the Higgs potential runs negative which would lead 
to a runaway instability at large field values. By intro- 
ducing Peccei-Quinn dynamics, we can potentially solve 
the strong CP problem, remove the unstable region, and 
obtain roughly the correct amount of dark matter due to 
a collection of statistical arguments that sets Fpq ^ E* . 
This is remarkably minimal, but does still leave questions 
regarding unification, baryogenesis, infiation, hierarchy 
problem, etc. It is conceivable that unification can still 
occur at higher energies by the introduction of new de- 
grees of freedom, that the physics of baryogenesis and in- 
flation is associated with such high scale physics [H, [3^ , 
and that the hierarchy problem has no dynamical expla- 
nation. Alternatively, the LHC or other experiments may 
discover new degrees of freedom at much lower energies, 
which would radically alter this picture. Currently all 
such issues remain largely unclear, requiring much more 
guidance from experiment and observation. 
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Appendix - Standard Model 2-Loop Beta Functions 

In this appendix we list the RG equations for the cou- 
plings A, yt,g' , g, gs at energies above the top mass rrit at 
2-loop order from Refs. In each case, we write 

dX/dt = /3\, etc, where t = In E/ii, and /i is the starting 
renormalization scale, taken to be mt ■ We also performed 
proper pole matching for couplings defined at the Z mass 
and running up to the top mass mt, but for brevity do 
not list those details here. 
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For the Higgs quartic coupling we have 



/3a = 



1 



(47r)2 
1 



6yt + I (2/ + {g' + g'^j') - {9g' + ig" I2yl) X 



(47r)4 
+A 



- (915.96 - 2895^5" - 559g2g'4 _ 3795'^) + SOyf ~ vt { ^ + i'2gl + 3A 



^9' + ^-^g'g" + ^-§9" + lOSg'X + - 312A^ 



79 



.4 , 21^2^,2 



19 



4^ 2 

For the top quark Yukawa coupling we have 
yt 



4 5'^ + A + ^5'^ + 80ff,2 - 144A 



19 



_9 2_ 17 2_ 2 9 2 
4 12 " 2 * 



(47r) 



23 4 3 2/2, 1187 ,4 



4 ^ ^5 .9 



216 



+ yr5s - 108.9^ + 



225 



131 



36g2 \y^ + 6 l^2yt - 2yi\ + X' 



16 " 16 

For the 3 gauge couplings gi — {g', g, gs} we have 



1 



{An} 



1 



(47r) 



where 



/4I 21 \ / 199/18 9/2 44/3 \ 

^= hr'-lT'-^ ' ^= 3/2 35/6 12 
V ^ ^ / \ 11/6 9/2 -26 / 



' 6 ' 2' 



By solving the set of 5 coupled differential equations, we obtain A as a function of energy or h. 



(16) 



(17) 



(18) 



(19) 



The wavefunction renormalization of the Higgs field is G{t) = exp(— '-f{t')dt'), where the anomalous dimension is 



1 



7 = 



(47r)2 
1 

(4^ 



95' ^ 35'2 „ 2' 
271 



431 



5 /9 



32 



16^ 



96 



2 V4 



17 



- i,9'9" - -^9" - ^ -y + ^9" + 85.^ + - 6A 



12' 



27 



(20) 
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